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ABSTRACT 


The  stiffiiess  and  strength  behavior  of  the  woven  fabric  composite  materials 
mainly  depend  on  the  properties  of  the  fiber  and  matrix  materials.  A  technique  was 
developed  to  design  and  analyze  woven  fabric  composites  based  on  the  micro-level 
properties  of  the  fiber  and  matrix  materials.  The  technique  is  based  on  the  bilateral 
relationships  among  the  fiber/matrix  materials,  the  unidirectional  composite  (strand),  the 
woven  fabric  layer  and  the  laminated  composite  structure  as  in  the  given  order. 
Simplified  and  efficient  analytical  models  were  developed  for  the  relationship  between 
any  subsequent  levels  of  woven  fabric  composite  materials.  The  technique  is  used  to 
predict  the  structural  level  stiffiiess  and  strength  in  terms  of  material  and  the  geometric 
configuration  of  the  woven  fabric  and  lamination.  Progressive  damage  or  failure  can  also 
be  simulated  at  the  fiber  and  matrix  level  by  using  this  technique. 
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1.  INTRODUCTION 


The  interest  in  composite  materials  has  been  growing  steadily  in  recent  years.  The 
use  of  composite  materials  in  the  fields  of  aeronautics  and  naval  construction  for  the 
fabrication  of  high  mechanical  performance  structures  is  increasing  day  by  day.  One  of 
the  most  commonly  used  composites  is  woven  fabric  composite,  in  which  strands  are 
formed  by  many  different  ways  of  weaving,  resulting  in  different  stiffiiess  and  strength. 
The  ease  of  handling  and  low  fabrication  cost  have  made  woven  fabric  composite 
materials  attractive.  However,  the  geometry  of  this  composite  class  is  complex,  and  the 
choice  of  architectures,  constituents  and  models  is  unlimited. 

Most  woven  composites  are  made  of  two  sets  of  interlaced  threads,  known  as  the 
warp  (vertical)  and  fill  (horizontal)  threads.  The  term  interlaced  is  used  to  indicate  where 
the  fill  and  warp  threads  cross  over  each  other.  The  warp  and  fill  threads  and  interlaced 
region  for  a  woven  composite  are  shown  in  Figure  1 . 

The  numerous  types  of  weaves  can  be  defined  by  the  repeating  patterns  in  both 
warp  and  fill  directions,  defined  by  the  n^g  and  The  number  of  n'^g  indicates  that  a 
warp  thread  is  interlaced  with  every  n^g  -  th  fill  thread  and  n^^g  means  that  a  fill  thread 
is  interlaced  with  every  n^'g  -th  warp  thread.  For  most  of  the  woven  composites  n^g  = 
n^g=ng.  Because  in  most  of  the  woven  fabrics,  only  one  type  of  fiber  material  is  used, 
material  parameters  like  thread  width  ratio  and  relative  fiber  volume  firaction  are  not 
used.  Woven  fabrics  with  n^4  and  with  interlaced  regions  that  are  not  connected  are 
known  as  satin  weaves.  As  defined  by  their  Ug  values,  the  fabrics  in  Figure  1  are  known 
as  plain  weave  (ng=2),  twill  weave  harness  satin  (ng=4). 
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Much  research  has  been  done  on  plain  weave  composites  particularly  on  the 
effective  stiffness  of  such  composites.  There  is  a  growing  need  to  develop  techniques  for 
the  analysis  of  woven  fabric  composite  structures  from  the  known  mechanical  properties 
of  constituent  phases.  In  order  to  predict  the  effective  material  properties  of  woven 
composites,  plate  lamination  theories  have  been  used  using  orthotrophic  yam  materials 
that  were  made  of  unidirectional  fibers  and  matrix.  [Ref.  1 -Ref. 8].  Less  attention  has  been 
given  to  the  strength  prediction  of  woven  composites. 

In  this  study,  a  new  multi-level,  micro/macro  modeling  for  a  woven  fabric 
composite  is  presented.  The  technique  can  be  used  to  predict  the  mayiTmim  stress  vmder  a 
certain  load  conditions  and  to  simulate  progressive  damage  of  a  woven  fabric  composite 
by  using  finite  element  formulation.  This  technique  has  a  similar  structure  to  the  one 
described  in  [Ref.  9].  However,  the  present  modeling  approach  is  much  simpler  and 
computationally  more  efficient. 
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II.  BACKGROUND 


A.  LITERATURE  SURVEY 

An  extensive  amount  of  research  has  been  conducted  to  gain  better  understanding 
of  woven  composite  materials.  This  is  important  because  the  improvement  of  woven 
fabric  composite  material  properties  depends  on  understanding  and  characterization  of 
the  basic  reinforcing  micro-mechanisms,  interaction  effects  between  the  various 
constituents  and  the  role  of  fiber  distribution  for  failure  progress. 

Ishikawa  and  Chou  [Ref  1]  developed  the  mosaic  model  that  was  based  on  the 
assumptions  of  constant  stress  and  constant  strain.  They  also  conducted  an  approximate 
analysis  by  talcing  into  account  the  fiber  imdulation  and  continuity.  Two  dimensional 
finite  element  analyses  support  their  mosaic  and  undulation  theory. 

Zhang  and  Harding  [Ref.  6]  tried  to  solve  the  problem  of  the  micro-mechanics 
analysis  of  the  elastic  constants  of  a  one-ply  plain  weave  composite  using  the  strain 
energy  equivalency  principle  with  the  aid  of  the  finite  element  method. 

Large  number  of  studies  was  conducted  to  determine  the  effective  material 
properties  fi-om  the  constituent  materials  and  their  geometry.  As  a  result  of  these  studies 
many  different  models  have  been  developed.  Most  of  these  studies  like  Moszisic  Model  of 
Ishikava  and  Chou  [Ref  2]  and  MESOTEX  of  Scida,  D.,  Aboura,  Z.,  Benzeggagh,  M.M., 
and  Boscherens,  E  [Ref  3]  focus  on  the  one-way  passage  from  the  constituent  material 
level  to  the  unidirectional  fiber  composite.  However,  the  other  passage  to  decompose 

m 

strains  and  stress  into  the  fiber  and  matrix  level  is  important  if  failure  or  damage  criteria 

are  to  be  applied  at  the  constituent  material  level.  Thus,  some  of  previous  works 

developed  the  bi-passage  models.[Ref  12-Ref  18].  Aboudi  [Ref  12]  and  Kwon  and  his 
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co-workers  [Ref.  13-Ref  15]  used  the  unit  cell  model  for  the  analytical  bi-directional 
module.  Aboudi’s  model  was  mathematically  elegant  and  Kwon’s  model  had  a  compact 
expression  for  computational  efficiency.  Pechold  and  Rahman  [Ref  16]  used  the 
strength  of  material  approach  with  equivalent  springs.  Their  model  was  very  simple  but 
neglected  the  interaction  of  Poisson’s  effects  of  the  fibers  and  the  matrix.  Hansen  and  his 
co-workers  [Ref  17-Ref  18]  used  another  approach  for  the  bi-directional  model. 
However  their  model  can  not  compute  the  effective  composite  material  properties  fi-om 
constituent  data  directly.  As  a  result,  they  used  an  independent  finite  element  analysis  at 
the  micro  level  to  determine  effective  composite  properties. 

Many  researchers  have  also  studied  the  relations  between  orthotropic  properties  of 
the  unidirectional  fiber-matrix  composite  and  the  effective  material  properties  of  the 
woven  fabric  composite.  In  order  to  predict  the  effective  material  properties  of  woven 
fabric  composite,  four  fimdamental  homogenization  methods  have  been  used  in  the  past: 
strain  energ>  balance,  plate  approximation,  direct  method  with  area  averaging  and 
asymptotic  expansion  homogenization. 

Strain  energ>'  balance  methods  exert  a  uniform  average  strain  while 
complementaiy  energy  methods  apply  a  uniform  average  stress.  The  asymptotic 
expansion  homogenization  approach  explicitly  employs  periodic  boundary  conditions  to 
simulate  the  condition  of  global  periodicity  of  the  microstructure.  Vandeurzen,  Ivens  and 
Verpoest  [Ref  9]  proposed  a  Combi  Cell  Method  (CCM)  for  a  general  weave  geometry 
which  is  capable  of  predicting  3D  stress  state  at  constituent  level  which  does  not  invoke 
isostrain  or  isostress  condition  throughout  the  unit  cell.  Gamich  and  Hansen  [Ref  17] 
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suggested  Multi-Continuum  Theory  (MCT).  This  bridged  the  structural  analysis  to  micro¬ 
mechanics. 
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III.  MULTILEVEL  TECHNIQUE 


A.  MODELLING  AND  FORMULATION 

The  proposed  multi-level  continuum  analysis  for  laminated,  woven-fabric 
composite  structures  is  sketched  in  Figure  2.  The  analysis  consists  of  three  modules.  The 
purpose  of  the  first  module  (called  Fiber-Strand  in  Figure  2)  is  to  relate  the  constituent 
materials  (i.e.  the  fiber  and  the  matrix  materials)  to  the  smeared  unidirectional  composite 
(strand)  material.  The  purpose  of  the  second  module  (called  Strand-Fabric)  is  to  connect 
the  unidirectional  composite  (strand)  material  to  the  woven  fabric  material.  The  third 
module  (called  Lamination)  is  designed  to  bridge  between  the  woven  fabric  material  and 
the  laminated  composite  structure  for  which  the  finite  element  analysis  is  performed. 

Each  module  has  bi-directional  passage  from  one  side  to  the  other.  In  other 
words,  each  module  computes  effective  material  properties  from  the  lower  level  material 
properties,  and  decomposes  the  upper  level  strains  into  the  lower  level  strains. 
Ultimately,  the  multi-level  analysis  determines  the  stresses  and  strains  at  the  fiber  and 
matrix  level  so  that  failure  or  damage  criteria  can  be  applied  to  the  constituent  materials. 

The  overall  procedure  of  the  multi-level  continuum  analysis  is  explained  below 
for  analyses  of  laminated  woven  fabric  composite  structures. 

1.  Compute  the  effective  material  properties  for  the  finite  element  structural  analysis 
using  the  three  modules  described  above,  starting  from  the  fiber  and  matrix 
properties. 

2.  For  a  load  increment,  compute  the  response  of  the  structure  such  as  deformation  and 
strains  from  the  finite  element  analysis. 
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3.  Apply  the  three  modules  from  top  to  bottom  to  compute  stresses  and  strains  at  the 
fiber  and  matrix  level. 

4.  Apply  failure  or  damage  criteria  to  check  potential  failure  or  damage  at  the 
constituent  material  level  such  as  fiber  breakage,  matrix  cracking,  and  interface 
failure. 

5.  Once  there  is  failure  or  damage,  degrade  the  material  properties  of  the  failed  or 
damaged  location  and  compute  the  effective  composite  properties  using  the  three 
modules  for  the  next  step  of  finite  element  structural  analysis.  If  there  is  no  failure, 
the  effective  properties  calculated  in  #1  can  be  used  for  the  next  finite  element 
analysis. 

6.  The  computation  continues  iteratively  until  the  composite  structure  fails  or  the 
designated  final  load  is  reached. 
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Figure  2.  Multi-level  Approach  to  Woven-Fabric  Composite  Structures 
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1.  Fiber-Strand  Module 

The  fiber-strand  module  relates  the  constituent  materials,  e.g.  the  fiber  and  matrix, 
to  the  smeared  unidirectional  composite  (strand).  Because  the  fiber  and  the  matrix  are  the 
base  units  for  composite  structures,  this  module  is  at  the  lowest  level  in  the  analysis.  The 
fiber-strand  module  has  two  functions;  the  first  function  is  to  compute  the  effective 
material  properties  of  the  unidirectional  composite  (strand)  using  the  fiber  and  matrix 
properties.  The  second  function  is  to  decompose  the  strains  and  stresses  at  the  strand 
level  to  those  at  the  fiber  and  matrix  level  so  that  progressive  damage  or  failure  can  be 
evaluated  at  the  fiber  and  matrix  level.  Thus,  the  Module  Fiber-Strand  establishes  the  bi¬ 
directional  relationship  between  the  constituent  materials  and  the  unidirectional 
composite  (strand). 

A  unit  cell  is  made  of  a  fiber  and  surrounding  matrix  material.  The  fiber 
represents  the  average  behavior  of  many  fibers  in  the  area  under  consideration,  and  is 
assumed  to  have  a  square  shape  of  cross-section.  Because  of  symmetry,  a  quarter  of  the 
unit  cell  is  considered  as  shown  in  Figure  3.  In  the  figure,  axis  1  is  the  fiber  direction, 
and  axes  2  and  3  are  normal  to  the  fiber  direction.  The  quarter  unit  cell  has  four  sub-cells 
and  the  size  of  the  sub-cell  depends  on  the  fiber  volume  fi:action. 

The  following  derivation  is  for  the  normal  components  of  stresses  and  strains. 
Stress  equilibrium  at  the  sub-cell  interfaces  is  expressed  in  terms  of  average  values  for 
each  sub-cell  as  below: 
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~  ^33  5 
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(3) 


where  subscripts  denote  stress  components  and  the  superscript  indicates  sub-cells,  as 
shown  in  Figure  3. 

Compatibility  of  deformation  is  assumed  for  strain  components  as  is  given  below: 


^11  ^11  ”^11’  ^11 

(4) 

44^ +  (1  -  4v^ )S22  =  4v^ £22  +  (1  “  4v^ )£22 

(5) 

4v^ £22 + (1  -  4v^  )£33 = 44^ £33 + (1  -  4i4  )s^2 

(6) 

4 +4  =  4 +4 

(7) 

^3  ^3  “  ^3  ^3 

(8) 

in  which  is  the  fiber  volmne  fi-action. 

The  unit  cell  (i.e.  strand)  stresses  and  strains  are  related  to  the  sub-cell  (i.e.  fiber 
and  matrix)  stresses  and  strains  as  below: 

=  v^cr|  +  4^  (1  -  4v^  )al  +  4v^  (1  -  4^  )<t|  +  (1  -  4v^  (jt  (9) 

=  v^e;  +4^ii-4^)s^.+  4i4(i  -  4^)s;  +  (i  -  44ff  (lo) 

These  equations  state  that  the  strand  stresses  and  strains  are  volumetric  averages  of  sub¬ 
cell  stresses  and  strains,  respectively. 

The  constitutive  equation  for  each  sub-cell  can  be  expressed  as: 

o;“  =  {ij,  k,l=\,2,ZwvA(x  =  a,b,c,d)  (11) 
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where  superscript  a  denotes  the  respective  sub-cell  of  fiber  or  matrix,  and  Efj^,  is  the 

material  property  tensor.  The  hygro-thermal  effects  are  not  included  in  Eq.  (11),  but  they 
can  be  easily  implemented  in  the  equation  as  shown  in  [Ref.  13]. 

Proper  algebraic  manipulations  of  the  previous  equations  finally  result  in: 

[£•']=  [K][£][S,]  (12) 

{e“}  =[«,]{«*}  (13) 

Where  [V]  is  the  matrix  containing  volume  fractions  of  each  sub-cell,  [£]  is  the  matrix 
consisting  of  material  properties  of  the  fiber  and  matrix,  and  [R^]  is  the  matrix  to 
decompose  the  strand  strains  into  the  fiber  and  matrix  strains  .  Equation  (12) 

computes  the  effective  material  property  matrix  for  the  strand  and  it  is  the  path 

from  the  constituent  level  to  the  strand  level. 

Conversely,  Eq.  (13)  is  the  path  from  the  strand  level  to  the  constituent  material 
level.  Therefore,  these  two  equations  establish  the  bi-directional  relationship  in  the  fiber- 
strand  module  and  the  fiber  and  matrix  stresses  are  computed  from  Eqs.  (1 1)  and  (13). 
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2. 


Strand-Fabric  Module 


The  strand-fabric  module  relates  material  properties  of  a  unidirectional  strand  to 
the  effective  material  properties  of  the  woven  fabric  composite.  Like  the  previous 
module,  this  module  has  two  functions:  computation  of  the  effective  properties  of  the 
woven  fabric  using  the  strand  materials  and  their  weaving  information,  and 
decomposition  of  the  woven  strains  and  stresses  into  the  strand  strains  and  stresses. 

Here,  the  strand-fabric  module  is  developed  for  a  plain  weave  composite.  The 
plain  weave  composite  model  is  presented  in  Figure  4.  The  model  has  13  sub-cells,  and 
most  of  those  sub-cells  have  fibers  along  the  axis-1  or  2  direction.  On  the  other  hand 
four  sub-cells  (#2,  4,  6,  and  8  in  Figure  4)  have  fibers  in  inclined  orientations,  and  sub¬ 
cell  #5  is  filled  with  the  matrix  material. 

First  of  all,  the  normal  stress/strain  components  are  discussed  below.  There  are 
three  normal  strains  for  each  sub-cell  with  the  total  of  39  strains  ,  n=\,  2, ...  13). 

There  are  three  normal  strains  at  the  unit-cell  level  ( ).  Those  strains  represent  the 
average  values  of  each  sub-cell  or  the  unit-cell,  respectively.  In  order  to  relate  them, 
equilibrium  and  compatibility  conditions  are  applied.  At  respective  interfaces  between 
any  two  neighboring  sub-cells,  the  normal  stresses  on  the  interface  plane  should  be  in 
equilibrium. 


,  {<■? = «)",  = {<r,  = «)^  (14) 

=  «)^  (O’  =«)’^  «)*  (O^  =«)'^  (15) 

(O’  =  (O^  =(0'“.  =(0%  (<^2T)^  =(0^  (16) 

(O'  =  (O'  (O'*  =  (0’',  (O'  =(0",  (17) 
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where  the  superscript  denotes  the  sub-cell  numbers  in  Figure  4  and  the  subscript  indicates 
the  stress  components.  In  addition,  some  strain  compatibility  conditions  are  assumed  for 


a  uniform  deformation. 

«)'+(«">"  =«)’  +  «)'“■  o») 

a,(.e‘T  +<>,(^.7)’  =<>.(0’  +4.^)'  +a,{0‘.  PO) 

+4,«D’  +«i(«,T)‘  “"iWD’  +4,(0*  +«,(«iT)'’.  (21) 

(£,r)’+(C)'’=(0’+(0'*.  (22) 

(4')'+(^")'*=(0’+(0'".  (23) 

+*,(4')'  +«2(«2?)’  =‘‘,(ssy  +42(«,t)’  +a,(^)\  (24) 

a,(4')=  +4,(4')’  +02(«S-)‘  =0,(4')*  +4,(4"')*  +“2(4')'*  (25) 

(4')’  +  (4')'’  =  (4')’+(4')",  (26) 

f,(4')'  +/:(0''’  =('.  +(2X0*.  (22) 

/,(4')’  +<,(4')"  =(/,  +»2)(0*,  (28) 

/,(4')>  +/2(4’)"  =(',  +(2X4')*.  (29) 

(£3^)' =(4')’.  (20) 

(4')*=(4')‘.  (21) 

',«)’  +(;(4')'*  =('l  +'2)(4')‘.  (22) 

(,(4)’  +  '2(4')'*  =  (',  +  '2)(4')*.  (23) 

(,(4')’  +(3(£3T)'*  =(',  +(2X4')*,  (34) 
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Here,  a, ,  b-, ,  and  are  the  dimensions  of  the  plain  weave  construction  as  shown  in 
Figure  4.  Finally,  the  umt-cell  strain  and  stress  are  assumed  to  be  the  volmne  average  of 
the  sub-cell  strains  and  stresses,  respectively.  That  is, 

(35) 

;j=l 

(36) 

n-\ 

Here,  F”  is  the  volume  fraction  of  the  n*  sub-cell.  The  constitutive  equation  for  each 
sub-cell  is  expressed  as: 

«)”  (37) 

Where  the  summation  sign  convention  is  applied  only  to  the  subscripts,  k  and  /.  For  each 
sub-cell,  should  be  determined  based  on  the  orientation  of  the  strand  inside  the 

sub-cell.  Algebraic  manipulation  of  the  previous  equations,  (i.e,  Eqs.  (14)  through  (37)) 
yields  the  following  relationship: 

=/((£, 5)".o„6,.(,)  (38) 

,a„b„t)  ■  (39) 

Equations  (38)  and  (39)  are  equivalent  to  Eqs.  (12)  and  (13)  for  the  fiber-strand 
module,  respectively.  Equation  (38)  computes  the  effective  material  properties  of  the 
woven  fabric  composite  based  on  the  material  and  geometric  properties  of  the  strands 
while.  Eq.  (39)  calculates  the  strains  at  the  strand  level  from  the  woven  fabric  strains. 
Once  the  strand  strains  are  computed,  strand  stresses  can  be  computed  from  the 
constituent  equations. 
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A  similar  derivation  can  be  made  for  shear  components.  For  example,  the  shear 
component  parallel  to  1-2  plane  can  be  derived  as  shown  below.  The  shear  stress 


equilibrium  at  the  sub-cell  interfaces  is  written  as 

^l(^lt)’+^2(^.TX'=(^.+^2X0-lTX. 

+?2«)''  =(^1  +^2X0% 

^,«0'’+^2«)''=(^+^2XO*, 

^,(0'+^2«)''=(^l+^2X0% 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 


Assumed  strain  compatibility  expressions  are: 


(4')’ =(«?)“> 

(O’ =(«.?)". 

«)’  =  («,T)’-  (51) 


Use  of  Eqs.  (40)  through  (51)  along  with  Eqs.  (35)  through  (37)  for  the  shear 
component  yields  the  shear  component  part  of  Eqs.  (38)  and  (39). 
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3. 


Lamination  Module 


The  lamination  module  eventually  bridges  the  effective  material  properties  of  a 
woven  fabric  lamina  and  the  finite  element  analysis  of  the  composite  structure. 
Depending  on  the  types  of  finite  elements  used,  the  connection  would  be  different  For 
laminated  composite  stmctures,  the  woven  lamina  properties  may  be  applied  to  the  finite 
element  analysis  program,  or  smeared  properties  of  several  laminates  may  be  pre¬ 
calculated  for  the  finite  element  program  using  a  lamination  theory.  In  this  study,  the 
shear  deformable  Mindlin-Reissner  plate  bending  theory  is  used.  The  theory  includes 
transverse  shear  deformation  and  calculates  the  smeared  material  properties  of  a  laminate. 

B.  FAILURE  CRITERIA  AT  THE  CONSTITUENT  LEVEL 

Failure  in  a  composite  stracture  is  described  in  terms  of  base  materials,  the  fiber, 
matrix  and  interface.  Thus,  failure  can  be  classified  as  fiber  failure,  matrix  failure  and 
interface  failure  at  the  constituent  material  level.  Since  fibers  are  long  and  slender, 
compressive  failure  of  the  fibers  is  quite  different  from  tensile  failure.  Thus,  four  failure 
are  considered:  fiber  failure  in  tension,  fiber  failure  in  compression,  matrix  failure,  and 
interface  failure. 

These  failure  modes  at  the  fiber  and  matrix  level  can  represent  various  modes  at 
the  structural  level  such  as  delamination,  fiber  splitting,  and  matrix  cracking.  Once  there 
is  a  failure,  the  corresponding  material  properties  are  degraded  accordingly.  Depending 
on  the  material  behavior  of  each  constituent  and  their  interface,  proper  failure  criteria 
need  to  be  selected.  For  example,  an  isotropic  failure  criterion  is  chosen  for  an  isotropic 
matrix  material. 
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VI.  NUMERICAL  SOLUTION  VERIFICATION 


First  of  all,  the  fiber-strand  module  and  the  strand-fabric  module  were  tested 
independently.  The  fiber-strand  module  was  validated  extensively  in  previous  works 
[Ref.  13- 15]  and  it  was  tested  again  for  a  carbon/epoxy  strand.  The  fiber  and  matrix 
materials  are  shown  in  Table  1.  The  effective  stiffness  and  strength  properties  of  the 
strand  made  of  the  fiber  and  matrix  materials  in  Table  1  and  with  a  fiber  volume  fraction 
of  0.7  are  compared  to  the  results  in  [Ref  19]  in  Table  2. 

Plain  weave  composites  were  studied  using  the  strand-fabric  module.  (Material 
properties  of  the  strands  of  the  woven-fabric  composites  are  tabulated  in  Table  3.)  The 
calculated  effective  material  properties  are  compared  to  other  results  in  Tables  4  and  5  for 
the  carbon/epoxy  and  e-glass/vinylester  woven  composites,  respectively.  The  present 
results  compared  well  to  other  experimental  data  or  analytical  results. 

Another  plain  weave  composite  made  of  the  strand  in  Table  2  was  also  studied  for 
its  tensile  strength.  The  predicted  strength  using  the  fiber-strand  and  strand-fabric 
modules  is  compared  to  the  test  result  in  [Ref  19].  As  shown  in  Table  6,  the  predicted 
strength  agrees  well  with  the  test  result 


Table  1:  Properties  of  Fiber  and  Matrix  Materials 


E2 

(GPa) 

Gi2 

(GPa) 

G23 

(GPa) 

Tn 

Tensile 

Strength 

(MPa) 

Shear 

Strength 

(MPa) 

Fiber 

221.0 

3.8 

13.8 

5.5 

20 

25 

3585 

4.4 

4.4 

1.6 

1.6 

0.34 

0.34 

159 

100 
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Table  2:  Properties  of  Strand  Material 


El 

(GPa) 

E2 

(GPa) 

Til 

ris 

Lx)ngitudinal 

Strength 

(MPa) 

Transverse 

Strength 

(MPa) 

Ref.  [19] 

151.0 

10.1 

0.24 

0.50 

2550 

152 

Present 

156.0 

10.2 

0.24 

0.54 

2543 

148 

Table  3 :  Material  Properties  of  Strands  and  Resin 


El  (GPa) 

E2(GPa) 

Gi2(GPa) 

G23  (GPa) 

/u 

T23 

E-GlassA^inylester 

57.5 

18.8 

7.44 

7.26 

0.25 

0.29 

Carbon/Epoxy 

134 

10.2 

5.52 

Ellfl 

EESi 

Epoxy 

\3A5 

3.45 

1.28 

1.28 

ES9 

Table  4:  Comi 

parison  of  P 

ain  Weave  Composite  Made  of  Carbon/Epoxy 

El  (GPa) 

E2(GPa) 

Gi2(GPa) 

G23(GPa) 

/l2 

Exp.  [Ref  20] 

55.5 

— 

4.93 

— 

0.0 

— 

Result  [Ref  5] 

56.1 

10.4 

5.08 

3.71 

Present  Result 
_ 1 

54.9 

10.2 

4.28 

3.47 

0.02 

0.47 

Table  5;  Com 

parison  of  Plain  Weave  Composite  Made  of  E-G] 

assWinylester 

E,  (GPa) 

E2(GPa) 

Gi2(GPa) 

G23(GPa) 

Yn 

r23 

Exp.  [Ref  3] 

24.8±1.1 

8.5  ±2.6 

6.5  ±0.8 

4.2  ±0.7 

Result  [Ref  3] 

25.3 

13.4 

5.19 

5.24 

0.12 

0.29 

Present  Result 

24.4 

12.4 

5.92 

4.52 

0.12 

0.28 

Table  6:  Plain  Weave,  Woven  Fabric  Composii 


Tensile  Strength  (MPa) 

753 

Ref  [19] 

750 

:e 
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V.  RESULTS  AND  DISCUSSION 


A.  PLATE  SUBJECTED  TO  IN  -PLANE  LOADING 

In  this  section,  square  shaped  woven  fabric  plates,  subjected  to  uniaxial  in-plane 

loading,  are  examined  for  fiber  failure  stress.  Because  of  symmetry,  one  quarter  of  the 
plates  was  modeled.  The  plates  are  plain  weave  composites.  The  fiber  and  matrix 
materials  are  listed  in  Table  1.  The  square  shaped  plates  are  meshed  with  100  linear 
triangular  elements.  The  fimte  element  mesh  can  be  seen  in  Figure  5.  The  longitudinal 
fiber  stresses  are  calculated  at  every  element. 


Symmetric  Boundary 

Figure  5.  Mesh  Generation  with  100  Linear  Triangular  Elements 
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1.  Effect  of  Fiber  Volume  Fraction  (V^)  on  The  Fiber  Failure  Stress 

Composite  material  properties  depend  upon  the  properties  of  fiber  and  matrix 
materials  of  which  they  are  made.  The  fiber  volume  fraction  is  the  percentage  of  the 
volume  of  composite  material  that  is  occupied  by  the  fiber.  Fibers  are  major  load- 
carrying  elements  in  a  structure  under  load.  As  a  result,  in  this  set  of  studies  the  effect  of 

fiber  volume  fi-action  on  the  fiber  stress  is  examined  for  the  in-plane  loaded  woven  fabric 
composites. 

Figure  6  shows  the  variation  of  the  maximum  fiber  stress  in  a  composite  as  a 
function  of  the  fiber  volume  fraction.  The  maximum  fiber  stress  is  normalized  with 
respect  to  the  value  when  the  fiber  volume  fraction  is  0.7.  Increasing  the  volume  fraction 
up  to  V^=0.7  increases  the  failure  stresses  linearly.  After  V^=0.7  increasing  results  in  a 
nonlinear  decrease  in  normalized  failure  stress. 


2.  Effect  of  Damage  Size  on  the  Fiber  Failure  Stress 

For  this  set  of  studies  the  elements  on  one  of  the  symmetric  edges  are  assumed  to 
be  damaged,  and  the  effect  of  damage  size  on  the  fiber  stress  is  examined.  The  damaged 
elements  have  the  reduction  of  material  properties  by  a  factor  of  100.  The 
fiber  stress  of  the  woven  composite  with  the  volume  fraction  V^=0.7,  and  no  damage 
condition  is  selected  for  normalization  of  the  computed  fiber  stresses.  Figure  7  shows 
the  fiber  stress  variation  with  along  with  the  damage  size  which  is  normalized  with 
respect  to  the  whole  length  of  the  symmetric  edge.  The  normalized  fiber  stress  shows  a 
non-linear  decrease  with  increasing  damaged  size. 
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B. 


PLATES  SUBJECTED  TO  BENDING  LOADS 


1.  Clamped  Plates 

In  this  section  of  the  study,  woven  fabric  plates  of  two  different  geometries 
(square,  rectangle)  with  the  same  surface  areas,  20x20,  10x40,  are  examined  for  the 
condition  of  clamped  edges,  under  bending  loads.  All  composite  plates  are  plain  weave 
composites.  The  fiber  and  matrix  materials  are  listed  in  Table  1  and  the  fiber  volume 
fraction  of  the  material  is  0.5.  Because  of  symmetry,  quarters  of  the  plates  are  modeled. 
The  quarter  plates  are  meshed  with  25  (5x5)  four-node  quadrilateral  elements.  The  mesh 
is  shown  in  Figure  8.  For  the  two  different  geometries,  plates  are  subjected  to 
concentrated  transverse  loads  at  the  center  and  distributed  pressure  loads,  respectively. 
The  longitudinal  fiber  stresses  are  computed  at  the  Gauss  integration  points  of  every 
element  due  to  the  fact  that  fibers  are  ultimate  load  carrying  members.  The  distribution  of 
fiber  stress  along  the  plates  can  be  seen  in  Figures  9-26.  The  fiber  stresses  are 
normalized  for  the  applied  unit  force. 

Plots  show  that  different  geometries  and  load  applications  result  in  different  fiber 
stress  values  and  stress  distributions.  For  the  square  plate  subjected  to  a  concentrated 
center  load,  the  maximum  fiber  stress  occurs  at  the  center,  and  the  minimum  fiber  stress 
occurs  at  the  comers  of  the  plate.  The  fiber  stress  decreases  fi-om  the  center  toward  the 
symmetric  boundaiy  until  the  2/3  distance  between  the  center  and  edge.  Then,  the  stress 
starts  to  increase  toward  the  symmetry  edge  where  the  fiber  stress  is  approximately  2/3  of 
the  maxirmim  fiber  Stress  at  the  center.  Due  to  the  symmetric  boundary  conditions  and 
equal  length  of  edges,  the  stress  distribution  is  symmetric  along  the  diagonal  direction  of 
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the  plate.  The  stress  distribution  for  the  concentrated  center  load  is  shown  in  Figures  9- 

10. 

Distributed  load  creates  the  maximum  fiber  stress  at  the  center  of  the  clamped 
boundaries.  The  fiber  stress  at  the  center  is  lower  than  the  stress  at  the  clamped  edge  but 
it  is  the  location  of  next  large  fiber  stress.  The  maximum  stress  caused  by  the 
concentrated  load  is  approximately  5  times  greater  than  the  maximum  stress  created  by 
the  distributed  load.  The  stress  distribution  for  the  distributed  load  can  be  seen  in  Figures 
11-12. 

For  the  rectangular  plate  subjected  to  a  concentrated  center  load,  the  maximum 
fiber  stress  occurs  at  the  center  again  while  minimum  stress  occurs  at  the  clamped  edges 
farther  away  from  the  center.  The  fiber  stress  decreases  from  the  center  toward  the 
clamped  edges  closer  to  the  center  imtil  the  midpoint  between  the  edge  and  center,  and 
then  it  increases  up  to  a  stress  value  at  the  edge,  which  is  very  close  to  maximum  stress  at 
the  center.  The  fiber  stress  is  negligible  on  about  the  half  of  the  plate  surface  near  the 
clamped  boundaries  farther  away  fi:om  the  center.  The  stress  distribution  for  the 
concentrated  center  load  is  shown  in  Figures  13-14. 

For  the  distributed  load,  the  maximum  fiber  stress  does  not  occur  at  the  center  but 
it  occurs  at  the  middle  of  long  clamped  edges.  The  fiber  stress  remains  almost  constant 
along  the  long  clamped  edges  until  the  half  of  the  distance.  Then,  the  stress  decreases 
almost  linearly  along  these  long  edges  and  it  becomes  zero  at  the  comers.  The  maximum 
fiber  stress  under  the  uniform  load  is  8  times  less  than  the  maximum  stress  under  the 
concentrated  load.  The  stress  distribution  for  the  distributed  load  can  be  seen  in  Figures 
15-16. 
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Clamped 


Fi^e  11.  Square  Plate,  Clamped  Edges,  Uniformly  Distributed  Load 
Fiber  Stress  Plot 
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Figure  12.  Square  Plate,  Clamped  Edges,  Uniformly  Distributed  Load 
Fiber  Stress  Contour  Plot 
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Figure  13.  Rectangular  Plate,  Clamped  Edges,  Concentrated  Load  at  The  Center,  Fiber 
Stress  Plot 


Figure  14.  Rectangular  Plate,  Clamped  Edges,  Concentrated  Load  at  The  Center,  Fiber 
Stress  Contour  Plot 
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Figure  IS.Rectangular  Plate,  Clamped  Edges,  Uniformly  Distributed  Load,  Fiber 
Stress  Plot 
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Figure  16.  Rectangular,  Clamped  Edges,  Uniformly  Distributed  Load,  Fiber 
Stress  Contour  Plot 


2.  Simply  Supported  Plates 

In  the  present  analysis,  the  plates  that  have  the  same  geometries,  material 
properties,  node  configurations  and  load  applications  as  the  plates  exammed  m  the 
previous  section  are  studied  for  simply  supported  edges. 

For  the  square  plate  subjected  to  a  center  concentrated  load  with  the  simply 
supported  edges,  the  maximum  fiber  stress  occurs  at  the  center  of  the  plate  where  the  load 
is  applied.  The  minimum  stress  occurs  at  the  midpoints  of  the  symmetric  edges.  Due  to 
the  symmetric  boundary  conditions  and  equal  length  of  the  edges,  the  stress  distribution 
shows  a  symmetric  behavior  with  respect  to  its  diagonal  of  the  plate.  The  stress 
distribution  for  the  concentrated  center  load  is  shown  in  Figures  17-18. 

For  the  distributed  load,  the  maximum  stress  occurs  along  a  region,  instead  of  just 
at  one  point.  This  region  is  just  like  a  canal  along  the  diagonal  of  the  plate.  The  stress  is 
distributed  symmetrically  with  respect  to  the  diagonal  of  the  plate  again.  The  stress 
distribution  for  the  distributed  load  is  shown  in  Figures  19-20. 

Rectangular  plates  with  simply  supported  edges  show  the  maximum  fiber  stress  at 
the  center  of  the  plates  for  both  distributed  and  concentrated  loads.  The  fiber  stress 
decreases  with  distance  from  the  center  in  all  directions,  and  finally  becomes  zero  at  the 
comers  and  midpoints  of  the  edges.  The  concentrated  load  results  in  faster  decrease  of 
the  fiber  stress  from  the  center  than  the  uniform  pressure  load. 

For  both  rectangular  and  square  plates,  the  concentrated  load  at  the  center  creates 
stress  values  that  are  five  times  greater  then  the  stress  values  created  by  distributed  loads. 
The  stress  distributions  are  shown  in  Figures  21-24. 
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Simply  Supported 
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Simply  Supported 


Figure  17.  Square  Plate,  Simply  Supported  Edges,  Concentrated  Load  at  The  Center, 
Fiber  Stress  Plot 


Figure  18.  Square  Plate,  Simply  Supported  Edges,  Concentrated  Load  at  The  Center, 
Fiber  Stress  Contour  Plot 
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Figure  19.  Square  Plate,  Simply  Supported  Edges,  Uniformly  Distributed  Load,  Fiber 
Stress  Plot 
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Simply  Supported 

Figure  20.  Square  Plate,  Simply  Supported  Edges,  Uniformly  Distributed  Load,  Fiber 
Stress  Contour  Plot 
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Figure  21.  Rectangular  Plate,  Simply  Supported  Edges,  Concentrated  Load  at  The  Center, 
Fiber  Stress  Plot 


Figure  22.  Rectangular  Plate,  Simply  Supported  Edges,  Concentrated  Load  at  The  Center, 
Fiber  Stress  Contour  Plot 
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Figure  23.  Rectangular  Plate,  Simply  Supported  Edges,  Uniformly  Distributed  Load, 
Fiber  Stress  Plot 
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Figure  24.  Rectangular  Plate,  Simply  Supported  Edges,  Uniformly  Distributed  Load 
Fiber  Stress  Contour  Plot 


C.  EFFECTS  OF  Ef,  AND  Gu  ON  THE  MAXIMUM  FIBER  STRESS 


In  this  section  the  effects  of  fiber  elastic  modulus  (Ef),  fiber  volume  firaction  of 
the  strand  (V  s)  3nd  the  undulation  angle  (Gu)  on  the  maximum  fiber  stress  are  examined. 
For  the  present  analyses,  square  plates  of  woven  fabric  composite  materials  with  clamped 
edges  and  the  concentrated  center  load  are  tested. 

As  the  elastic  modulus  of  the  fiber  is  increased  gradually,  maximum  fiber  stresses 
are  computed.  The  result  is  plotted  in  Figure  25.  The  graph  is  normalized  with  respect  to 
the  elastic  modulus  of  the  fiber  in  Table  1,  and  the  maximmn  stress  value  in  Figure  9, 
respectively.  As  the  fiber  modulus  increases  by  10  times,  the  maximum  fiber  stress 
increases  more  than  40  percent. 

Because  the  fibers  are  major  load  canying  members,  higher  fiber  modulus  makes 
fibers  carry  more  stresses  with  decreasing  stresses  on  the  matrix  material.  The  mavimnm 
fiber  stress  increases  nonlinearly  with  the  increasing  fiber  modulus.  As  the  fiber  modulus 
to  the  matrix  modulus  ratio  increases,  the  effect  of  fiber  modulus  is  smaller  on  the 
maximum  fiber  stress. 

Figure  26  shows  the  fiber  stress  variation  as  a  function  of  the  fiber  volume 
fi'action  of  the  strand.  The  maximum  fiber  stress  decreases  with  a  smooth  nonlinearity 
with  increasing  fiber  volume  fraction. 

Figure  27  shows  the  maximum  fiber  stress  variation  in  terms  of  the  undulation 
angle  of  the  plain  weave  composite.  Undulation  angle  is  the  angle  that  a  yam  makes  with 
the  horizontal  axis  at  the  interlaced  regions.  The  higher  the  undulation  angle  is,  the  higher 
the  maximum  fiber  stress  becomes.  The  relationship  is  approximately  linear. 
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D.  FAILURE  PROGRESS 

In  this  section  of  the  study,  the  sqiiare  plates  studied  in  the  previous  sections  are 
examined  for  the  failure  progress.  Most  of  the  plates  show  that  when  the  failure  initiates 
at  the  point  of  plate  bearing  the  greatest  fiber  stress,  such  failure  progresses 
catastrophically  with  load  transfer  between  the  failed  zones  and  adjacent  zones.  Once  the 
adjacent  zones  carry  an  extra  load  due  to  the  failed  zone,  the  neighboring  zones  start  to 
fail. 

For  the  square  plate  with  simply  supported  edges  and  uniformly  distributed  load, 
as  shown  in  Figures  28  and  29,  the  maximum  fiber  stress  is  experienced  along  the 
diagonal  direction  between  the  center  and  the  comer.  When  this  fiber  stress  value  reaches 
the  failure  strength,  it  will  cause  the  stracture  to  fail  along  the  diagonal. 
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Simply  Supported  Simply  Supported 

Figure  28.  Maximum  Fiber  Stress  Occurs  at  the  Center  and  Failure  Initiates  at  That  Point 
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Simply  Supported  Simply  Supported 

Figure  29.  Failure  at  The  Center  Gets  Bigger  Catastrophically  and  The  Whole  Structure 
Fails 
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E.  SHELL  STRUCTURES  SUBJECTED  TO  BENDING  LOADS 


In  this  section,  cylinder-shaped  woven  fabric  shell  structures,  subjected  to 
concentrated  loads,  are  examined  for  fiber  stress  distribution  along  the  structure.  Fiber 
and  matrix  materials  are  the  same  as  in  the  previous  sections,  and  fiber  volume  firaction  is 
0.5.  For  this  set  of  studies,  one  edge  of  the  cylinder  is  clamped  and  one  of  the  comers  on 
this  edge  has  a  notch.  The  notch  is  of  square  shape  and  its  size  is  1/10  of  the  cylinder 
edge.  The  stmcture  is  meshed  with  400  (20X20)  four-nodded  quadrilateral  elements.  The 
concentrated  load  is  applied  at  3  different  points.  The  mesh  generation  and  load 
application  points  are  shown  in  Figure  30.  For  three  different  points  of  load  application, 
the  fiber  stress  distributions  along  the  stmcture  can  be  seen  in  Figures  31-33. 

The  load  applied  at  the  comer  just  above  the  notch  (point  PI)  creates  the 
maximum  fiber  stress  at  the  comer  of  the  notch.  The  fiber  stress  decreases  with  distance 
until  reaching  the  3/5  of  the  clamped  edge  and  increases  again  towards  the  other  comer  of 
the  clamped  edge.  The  fiber  stress  at  load  application  point  is  1/5  of  the  mavimiim  fiber 
stress  at  the  notch.  The  stress  distribution  can  be  seen  on  Figure  31. 
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Figure  30.  Cylinder  Shaped  Shell  Structure  Meshed  with  400(20x20)  Equilateral 
Elements 
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Figure  31.  Fiber  Stress  Distribution  Along  The  Cylinder-Shaped  Structure,  Bending  Load 
at  Point  PI 
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The  load  at  the  center  of  upper  edge  (point  P2)  again  creates  the  maximum  fiber 
stress  at  the  comer  of  the  notch.  The  fiber  stress  at  the  load  application  point  is  1/2  of  the 
mavirmim  fiber  stress.  The  fiber  stress  decreases  with  distance  fi:om  the  notch  comer,  and 
the  stress  is  distributed  symmetrically  around  the  load  application  point.  The  stress 
distribution  can  be  seen  in  Figure  32. 

The  load  applied  at  the  comer  furthest  from  the  notch  (point  P3)  creates  the 
mavirmim  stress  at  the  closest  point  to  the  load  along  the  clamped  edge.  The  fiber  stress 
at  the  notch  comer  is  1/2  of  the  maximum  fiber  stress  and  the  fiber  stress  at  the  load 
application  point  is  1/5  of  the  maximum  fiber  stress.  The  fiber  stress  distribution  is  shown 
on  Figme  33. 

The  mavimiim  stress  value  and  location  change  according  to  the  load  application 
points.  As  the  loads  get  closer  to  the  notch  comer,  the  maximum  fiber  stress  increases 
around  the  notch. 
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Figure  32.  Fiber  Stress  Distribution  Along  The  Cylinder-Shaped  Structure,  Bending  Load 
at  Point  P2 
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Figure  33.  Fiber  Stress  Distribution  Along  The  Cylinder  Shaped  Structure,  Bending  Load 
at  Point  P3 
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VI.  CONCLUSION  AND  RECOMMENDATIONS 


A  multi-level  technique  for  woven  composite  materials  is  presented  in  this  study. 
The  technique  can  be  used  for  simulation  of  progressive  failure  in  woven  fabric 
composite  to  predict  residual  stiffness  and  strength  along  with  the  progressive  failure. 
The  results  of  effective  stiffiiess  and  strength  computations  for  woven  composite 
materials  with  multi-level  technique  are  in  good  agreement  with  data  obtained  from 
experiments  and  other  finite  element  analyses. 

Woven  fabric  composite  plates  subjected  to  in-plane  loads  and  bending  loads  and 
woven  fabric  shell  structures  subjected  to  bending  loads  are  studied  parametrically.  The 
longitudinal  fiber  stresses  are  computed  using  the  multi-level  technique,  because  those 
fibers  are  major  load  carrying  elements.  The  important  results  of  parametric  studies  for 
the  variation  of  fiber  stresses  are  summarized  in  the  following.  These  results  demonstrate 
the  usefulness  of  the  proposed  multi-level  approach  for  computational  modeling  and 
simulation  of  woven  fabric  composite  structures. 

1.  Plates  with  different  geometrical  shapes  and  load  applications  result  in  different 
mavirmim  fiber  stress  values  and  fiber  stress  distributions.  Concentrated  loads  always 
create  greater  tnaYimum  fiber  stress  values  than  the  maximum  fiber  stress  values  created 
by  distributed  loads. 

2.  An  increase  in  elastic  modulus  and  undulation  angle  increases  the  maximum  fiber 
stress.  However  an  increase  in  fiber  volume  fraction  decreases  the  maximum  fiber  stress. 

3.  For  the  plate,  when  failure  is  initiated  at  a  certain  part  of  the  plate,  bearing  the 
greatest  fiber  stress,  failure  progress  catastrophically  with  load  transfer  between  the  failed 
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zones  and  adjacent  zones.  Once  the  adjacent  zones  cany  extra  load  due  to  the  failed  zone, 
the  neighboring  zones  begin  to  fail  again 

4.  For  the  shell  structures,  which  have  notches,  the  maximum  stress  value  and 

location  change  according  to  the  distance  between  the  load  application  point  and  the 
notch. 

This  study  can  be  extended  to  model  other  t5^es  of  woven  fabric  composites  such 
as  twill  weave  and  four  harness  satin  woven  composites.  Different  plate  and  shell 
structures  may  be  analyzed  with  different  boundary  conditions.  Some  experimental  work 
is  recommended  to  evaluate  the  analytical  study  conducted  here 
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